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Abstract 

Using Stein's method and the Malliavin calculus of variations, we derive explicit estimates for the 
Gamma approximation of functionals of a Poisson measure. In particular, conditions are presented 
under which the distribution of a sequence of multiple Wiener-Ito stochastic integrals with respect 
to a compensated Poisson measure converges to a Gamma distribution. As an illustration, we 
present a quantitative version and a non-central extension of a classical theorem by de Jong in 
the case of degenerate [/-statistics of order two. Several multidimensional extensions, in particular 
CNI ' allowing for mixed or hybrid limit theorems, are also provided. 
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1 Introduction 

The use of the Malliavin calculus of variations in order to deduce limit theorems for non-linear func- 
tionals of random measures has recently become a relevant direction of research, one reason for that 
being the many successful applications in geometric probability or stochastic geometry. Apart from 
a few exceptions, most contributions to this topic fall into the two categories of normal and Poisson 
approximations; see [21 [TOj, HU 031 E21 [Ml EH [38] for distinguished examples of the former class, mostly 
qs, • based on the use of the Stein's method (cf. [TS]); see [2l[2TJ[3l] for references based on the combination 
of Malliavin calculus and of the Chen-Stein method for Poisson approximations. We also refer to [7] 
for recent extensions to general absolutely continuous distributions having support equal to the real 
line. 



The aim of the present paper is to provide the first array of results concerning limit theorems on 
the Poisson space, where the limit distribution is absolutely continuous and has support contained 
in a proper subset of R. More precisely, we are interested in probabilistic approximations where the 
limiting random variable has a centred Gamma distribution T v with parameter v > 0. We say that a 

random variable G{y) has distribution T u if G{y) = 2F(v/2) — v, where F{y/2) has a usual Gamma 

distribution with mean and variance both equal to v/2 (here and throughout = stands for equality in 
distribution). If v > 1 is an integer, then T u reduces to the centred ^-distribution with v degrees of 
freedom. We remark that the support of T u is given by the half-line [— u, +oo), and that the first four 
moments of T u are 0, 2v, 8u and 12z/ 2 + A8is, respectively. We will often meet these expressions in the 
discussion to follow. 

Our main contribution is the general estimate stated in Theorem 12.14 which involves Malliavin 
operators and is obtained by means of Stein's method, allowing one to measure the distance between 
the law of a given Poisson functional and T u . This estimate is applied to deduce explicit sufficient 
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conditions for Gamma limit theorems involving sequences of multiple Wiener-Ito stochastic integrals. 
Our analysis is significantly inspired by [15j [16] , where the problem addressed in the present paper 
was first studied in the framework of non-linear functionals of general Gaussian fields. However, due 
to the combinatorial complications one has to face when dealing with point measures, our paper con- 
tains a number of new subtle computations related to the explicit estimation of Malliavin operators 
on configuration spaces. One specific problem we will have to deal with is that the solution of the 
Stein's equation associated with the law of G{v) is not differentiable at x = —v. Thus, in order to 
obtain bounds that are well-suited for our applications (which may involve random variables possibly 
taking values in (—00, —v)), we will have to combine techniques recently introduced by Schulte [33J 
with classical isometric formulae borrowed from the standard reference [27]; see Proposition 12.31 be- 
low. One should note that, in view of the exact chain rules that are available on a Gaussian space, 
the non-differentiability of the Stein solution in one point is immaterial when studying the Gamma 
approximation of smooth functionals of a Gaussian field; see again [15|. 116]. 

As an illustration, we will include some applications to non-central limit theorems for sequences of 
degenerate (in the sense of Hoeffding) [/-statistics. Our findings generalize several classic result in the 
field; cf. [H |8]. In particular, we derive a quantitative and a non-central version of a famous theorem 
by P. de Jong [U [5] . Our analysis also contains a quantitative version of a non-central result recently 
discussed by Reitzner and Schulte [301 Section 5.1]. 

Finally, to demonstrate the flexibility and scope of our approach, we will show that our analysis 
can naturally be extended to a multidimensional framework. We will not only obtain multidimensional 
Gamma limit theorems, but also mixed ox hybrid results, where the multidimensional limit distribution 
is composed both of Gamma and of normal or Poisson components. This kind of limit theorems heavily 
relies on our use of Malliavin operators. We are not aware of any other available technique allowing 
one to deduce general mixed limit results, such as the ones deduced in the present paper. We shall see 
that our findings are a refinement of the 'Portmanteau inequalities' recently obtained by Bourguin and 
Peccati in [2]. In this respect, we stress that our results will implicitly yield a collection of sufficient 
conditions in order to have that two sequences of Poisson functionals are asymptotically independent. 
This provides a new contribution to the difficult and mostly open problem of characterizing the 
asymptotic and non-asymptotic independence of functionals of a Poisson measure; see e.g. |28t 131]. 

The remainder of the paper is organized as follows. In Section [2] we present our results in full 
generality. Some background material is collected in Section (3[ whereas the final Section [3] contains 
detailed proofs, as well as some ancillary technical results. 

2 Presentation of the results 

We will now present an overview of the main findings of the paper. To enhance the readability of 
our text, we have gathered together in Section definitions, notation and relevant results from the 
literature. 

2.1 General limit theorems 

Every random object considered below is defined on a suitable probability space (tljJ 7 , P). The 
approximation results obtained in the present paper deal with (real-valued) functionals of a Poisson 
measure rj on some Polish space (Z, 3f) having non-atomic and cr-finite control /U; see Section [3]- (I). We 
will assume that these functionals are square- integrable random variables. To measure the distance 
between the distribution of a functional F of ij and that of a centred Gamma random variable G(u), 
we shall use the (pseudo-) metric ^3, which is defined as follows: for every pair of square-integrable 
random variables X, Y, we put 

d 3 (X,Y) = sup \E[h(X)]-E[h(Y)}\, 

hen 3 
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where rl 3 := {h E C 3 : [|/i^[|oo < 1) i € {1,2,3}} (with faw the derivative of order j of /i), and where 
C 3 is the space of thrice differentiable functions on M having bounded derivatives. We notice that the 
topology induced by d 3 is stronger than the topology induced by convergence in distribution, which 
implies that if d 3 (F n , G{yy\ — >• 0, as n — > oo, for some sequence of functionals F n , then the distribution 
of F n converges to T v . By a slight abuse of notation, and to stress the role of the underlying Gamma 
distribution, we shall often write d 3 (F, T v ) instead of d 3 \F,G(v)\. 

For q > 1, we write L 2 (/i q ) to indicate the Hilbert space of Borel-measurable functionals on Z q that 
are square-integrable with respect to \i q . We also use the following special notation: L 2 ^ 1 ) = L 2 (/i), 
and L 2 ym (n q ) is the subspace of L 2 (/x 9 ) composed of those functions that are fj, q -a.e. symmetric; see 
Section [HJ- (II). Moreover, in order to simplify the notation, we use the convention that || • || and ( • , • ) 
stand for the norm and the scalar product in some space L 2 (^, q ) whose order q will always be clear 
from the context. 

Our first result is a quantitative estimate for d 3 (i 7 , T u ^j in terms of the Malliavin operators D and 
L _1 , that is, the derivative operator and the pseudo-inverse of the Ornstein-Uhlenbeck generator. We 
recall that the derivatives DF and DL~ 1 F are random elements with values in the Hilbert space 
L 2 (fi); see Section d(V). 

Theorem 2.1 (General Gamma bounds). Let F be a centred and square-integrable functional of 
the Poisson measure ij, and assume that F is in the domain of the derivative operator D. Then, 

(2.1) d 3 (F,T u ) < c 1 A 1 (F) + c 2 A 2 + 2c 1 A 3 (F) 

:= ciE \2(F + v) + - (DF, -DL,- l F) \ + c 2 J ^D.F^D.L- 1 F\] fi(dz) 

+ 2c 1 ^E[(Z) 2 l {F> _, } )( J D z F)|^L- 1 F|] fi(dz), 

with constants c\ and c 2 given by 

a = max(l,l/z/ + 2/z/ 2 ) and c 2 = max(2/3, 2/(3^) - 3/V 2 + 4/i/ 3 ). 
If in addition E[(DF, -DL~ 1 F)\F] >0 (a.s.-P), then 



A^F) < A[(F) := (2(F + v)- (DF, -DL^F)) 2 ] , 

and consequently 

(2.2) d 3 (F,T u ) < Cl A[(F) + c 2 A 2 + 2 Cl A 3 (F) 

Remark 2.2. (i) In (|2,ip . we implicitly used a 'trajectorial' definition of the random function 
z h-> D z l{ F> _ u y, that is, we put D z l{ F> _ u y = l^F + £, z p > _ u ^ — l^p > _ u y, without necessarily 
assuming that 



E 



JjD z l {F> _ u} ) 2 n(dz) < oo 



(note that this last relation is equivalent to the fact that l{p>-v} belongs to the set domD, as 
defined in Section 3-(V); see Lemma 133]) . It is easily checked that 

{D z l{F > -v})(D z F) = (l{ F <_ v<F+DzF } + 1{F+DzF<-v<F})\D z F\, 

in such a way that A 3 (F) > 0. An effective bound on A 3 (F), in the case where n is a finite 
measure and F is a multiple Wiener- ltd integral, is presented in Proposition [ 



(ii) As first done in [22], we shall often control the quantity A 2 (F) appearing in (|2.ip by using the 
relation 

(2.3) A 2 (F) < A 4 (F) x A 5 (F) := ^jf E[|I^F| 4 ] /i(dz)) ' x ^jf E[\D Z L' X F\ 2 } p(dz)) ' . 
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We also note that, if {F n : n > 1} is a sequence of random variables with bounded variances 
living in a fixed sum of Wiener chaoses, then the numerical sequence n h-> A^{F n ) is necessarily 
bounded. 

(iii) Theorem 12.11 should be compared with the following bound from |16( Theorem 3.11]. Let F be 
a centered functional of a Gaussian measure on Z with control fi, and assume that F is in the 
domain of the Malliavin derivative D (see |18|. Chapter 2] for relevant definitions), then there 
exists a constant K such that, for some adequate distance d, 

d(F,T v ) < K x E \2{F + - {DF, -DL~ 1 F)\ . 

The presence of the additional term 

c 2 f E[\D z F\ 2 \D z L~ 1 F\}^dz) + 2c x f E[{D z l {F> _ u} ){D z F)\D z L-- x F\] fi(dz) 
Jz J z 

in (|2.1|) or (|2.2p is due to the characterization of the Malliavin derivative on the Poisson space 
as a difference operator as well as to the non-differentiability at — v of the solution of the 
Stein-equation characterizing T u ; see Section [3]- (V). As proved in [TBJ Proposition 3.9], on the 
Gaussian- Wiener space the condition E [(DF, -DL^F^F] > ( a.s.-P) is satisfied for every F 
in the domain of D. 

(iv) Other relevant one-dimensional bounds for probabilistic approximations involving Malliavin op- 
erators on the Poisson space are proved in [22], dealing with normal approximations, [21], dealing 
with the Poisson approximation of integer- valued random variables and [7j, focusing on abso- 
lutely continuous distributions whose support is given by the real line. See [2j [23] for several 
multidimensional extensions. 

As announced, we conclude the present section with a useful bound on the quantity A^(F), in the 
case where F = I q (f) equals a multiple Wiener-Ito integral and the control measure [i is finite. At the 
cost of a heavier notation, our techniques could suitably be modified in order to deal with the case of 
a random variable F having a finite chaotic expansion. 

Proposition 2.3. Let the control measure \x be finite, and consider F = I q (f), where q > 2 and 
f ^ Fg ym (fi q ). We assume that (i) E j z (D z F) A fi(dz) < oo, that (ii) the random function 

Z 3 z ^ D Z F\D Z F\ := v(z) 

is such that v(z) £ domD for fi(dz)-almost every z, and satisfies 

Elf (D Z2 v(*i)) 2 »{dzi)v(dz 2 ) < oo. 
Jz Jz 

Then, defining A^(F) as in \2. one has the bound 



2^2 

(2.4) 



A 3 (F) < W;E^(^F)V(d^) + W;E^^(Z) 22 I) 2l F) 2 p 2l F)V(d2i)/i(dz 2 ) 
+ JeJ^J (D Z2 D zl F)^{d Zl )fx(dz 2 ). 



Remark 2.4. (i) Another way of controlling the term A^(F), whenever F has a finite chaotic ex- 
pansion, is discussed in [33]. One should note that, albeit our proof of Proposition 12.31 also starts 
with an integration by parts formula, our strategy for controlling the term ^3 (F) is significantly 
different. Indeed, our approach is based on isometric formulae for divergence operators, whereas 
|33j uses a direct estimation consisting in controlling \DF\ by a random function having a finite 
chaotic expansion. When applied to our framework in the case q > 2, the technique used in 
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|33j leads to expressions involving contractions of the absolute value of the kernel /, therefore 
producing bounds that are systematically larger than ours. When applied to the case q = 2, the 
strategy adopted in [33] leads to slower rates of convergence, but allows in principle to dispense 
with the assumption that the underlying cont rol measure has finite mass. Since all our applica- 
tions concern sequences of control measures having a finite mass, and for the sake of conciseness, 
we will omit a formal discussion of this fact. 

(ii) From the standpoint of geometric applications, focusing on Poisson measures having a finite 
control is barely a restriction. Indeed, the kind of geometric limit theorems we are interested 
in typically involve either functionals of a Poisson measure having a finite control, whose total 
mass asymptotically explodes (like the ones we consider in the applications developed later in 
the paper) , or functionals of the restriction of a Poisson measure to a finite window with growing 
volume; see e.g. [21 El [TOj, [Til E2J Ell EH EI] for a recent collection of distinguished examples. 

2.2 Simplified estimates for supports contained in a half-line 

The applications we are interested in require that the we consider random variables possibly taking 
values in the half-line (— oo, — u), in such a way that the rather unusual term A 3 (F) cannot be dispensed 
with. However, if one is only interested in measuring the distance between and the law of a random 
variable with support in [— u, +oo), then the statement of Theorem 12.11 can be significantly simplified, 
since in this case the term A 3 (F) disappears. In particular, whenever the law of F satisfies these 
requirements, the finiteness of the measure \i does not play any role. This point is made clear in the 
next statement whose easy proof is left to the reader. 

Proposition 2.5. Let F be a centered square-integrable functional of the random measure ij. Assume 
that the law of F has support in [— u, +oo) and that F is in the domain of the derivative operator D. 
Then, the bound l2l\) holds with A 3 (F) = 0. // moreover E[(DF, -DL~ 1 F)\F] > (a.s.-F), then the 
estimate l2Jfy holds with A 3 (F) = 0. 

2.3 General results for sequences of multiple integrals 

We now focus on the following setup. Let (Z, be a fixed Polish space as above, and {r] n : n > 1} 
be a sequence of Poisson random measures on (Z, 2f), such that, for each n, the non-atomic control 
measure [i n of rj n is finite. In view of applications, we allow that /x„(iJ) — > oo, as n — > oo. For a given 
even integer q > 2, we consider a sequence {I q (f n ) : n > 1} of multiple Wiener-Ito stochastic integrals 
with the following characteristics: (a) {/„ : n > 1} C L 2 ym (/in) is composed of kernels satisfying 
the technical assumptions stated in Section [31- (VIII) below, and (b) for every n > 1, the integral 
Iq(fn) is realized with respect to the compensated Poisson measure fj n = rj n — fi n . The next theorem 
characterizes the convergence of the distribution of I q (f n ), as n — > oo, to the limit law T u . The set 
of analytic conditions appearing below is expressed in terms of (possibly symmetrized) contraction 
kernels, whose definition is provided in Section [3]- (VI). Observe in particular that f n -k q f n = f%. 

Theorem 2.6 (Gamma limits in the Poisson- Wiener chaos). Let the above assumptions and 
notation prevail (in particular, [i n is a finite measure for every n), let q > 2 be an even integer and 
let {f n : n > 1} C L\ vm (nn) be such that lim (?!||/ n || 2 = 2^, and suppose that the technical conditions 
of Section^ (VLII) are satisfied. Assume in addition that 

(2.5) lim || f n -k l r f n \\ =0 and lim \\f n l?J/tf n - c q f n \\ = with c q = 4 2 

for all pairs (r, £) such that either r = q and £ = 0, or r € {1, . . . , q}, I E {1, . . . , min(r, q — 1)} and 
r and i are not equal to q/2 at the same time. Then, the distribution of I q (f n ) converges to T u as 
n — > oo. Moreover, for some positive finite constant K independent of n, 

d 3 (I q (f n ),T u ) < dA^fn)) + C 2 A 4 (L q (f n )) X A 5 (I q (f n )) + 2 Cl A 3 (I q (f n )) 

(2.6) < K x maxjl^HMI 2 - 2i/|; ||/„ *£ f n \\; \\f n 4 /n|| 1/2 ; H/n^/n " c q f n \\} -> 0, 
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where we have used the notation introduced in \2.1\) - l2lfy . and the maximum is taken over all p = 
1, . . . , q— 1 such that p / q/2 and all (r, £) such that r / £ and either r = q and I = ; or r E {1, . . . , q} 
and i 6 {1, . . . , min(r, q — 1)}. 

Example 2.7. (i) Assume q = 2. Then, C2 = 1 and the maximum in (|2,6p is taken over the 
following four quantities: 

\2\\f n \\ 2 - 2v\, ||/„^/ n || 1/2 , \\fn4fn\\ 1/2 , ||/n*i/n-/n||- 

(ii) Assume q = 4. Then, C4 = 1/18 and the maximum in (|2.6p is taken over the following ten 
quantities: 

|2||/„|| 2 -2Z,|, \\fn4fn\\, ||/n^/n|| 1/2 , \\fn 4 fn\\ 1/2 , \\fn 4 fn\\ l/ \ ||/n ^ /n|| 1/2 , 
||/n^/n|| 1/2 , ||/n^/n|| 1/2 , || fn 4 fn II ^ , \\fn *\fn ~ IS^/nll , 

where we have used the fact that ||/ n *i /n|| = \\fn *| /nil- 
Remark 2.8. (i) Under the assumptions in the statement, one has that the sequence 

1/2 



MWn)) ■= (^J z nD z L- l I q (f n )\ 2 ]ll n {&z) 



A 5 (I q (f n )) 2 = (q - l)\\\f n \\ 2 ^ — >0 as n -> oo. 



is such that 



It follows that our inequality (|2.6p not only provides an analytic bound in the distance ds, but 
also ensures that the three numerical sequences {Ai(I q (f n )) : n > 1}, {A^{I q {f ri )) : n > 1} and 
{A^(I q (f n )) : n > 1} (all related to Malliavin operators) converge to zero. This fact is crucial 
when dealing with the multidimensional results discussed in Section 12,61 An analogous remark 
applies to Proposition 12.91 and Theorem 12.131 below. 

(ii) Similar conditions (only involving contractions of the type *£, with r = 1, . . . ,q — 1) in the case 
of multiple integrals with respect to a Gaussian measure can be found in |15|. Theorem 1.2]. 
Non-central results of a similar flavor, in the context of free probability and multiple integrals 
with respect to a free Brownian motion, are proved in [17) . 

(iii) We were able to deduce meaningful conditions for Gamma approximations only in the case of 
an even integer q > 2. However, unlike in the Gaussian case (see |15l Remark 1.3]), in a Poisson 
framework one cannot exclude a priori the existence of a sequence of multiple integrals of odd 
order converging to a limiting Gamma distribution. We prefer to consider this issue as a separate 
problem, and keep it as an open direction for future research. 

(iv) In the estimate (j2.6|) . and in contrast to the main bounds on normal approximations proved in 
[2"2"j . norms of the type \\f n 4 f n \\, r ^ £, appear under a square root. This phenomenon seems 
unavoidable, and it is directly related to the presence of cross terms arising from the specific 
form of the Stein equation associated with the Gamma distribution. 

The following statement shows that condition (|2,5p might take a particularly attractive form in the 
case of double Poisson integrals. This will be used in order to prove the results presented in Section 
2.41 dealing with the Gamma approximation of degenerate [/-statistics. 



Proposition 2.9 (Three moments suffice for Gamma approximations). Let the control meas- 
ures {fi n : n > 1} be finite, let q = 2 and let {f n : n > 1} C L 2 (/i 2 ) be such that lim E[/ 2 (/ n )] = 

lim 2||/ n || 2 = 2u, and such that the technical conditions of Section^ (VIII) are satisfied. Assume in 
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addition that f z f^dfi^ — > and that Efi^C/n)] < 00 f or every n. Then, condition \2. 5\) is verified if 
and only if 



(2.7) 



nit(fn)} - 12E[/f (/„)] — >■ 12i/ 2 - 48^ as 



n 



00. 



In particular, if the sequence is uniformly integrable, then \2. 5\) and (2.1) are both necessary and 
sufficient in order to have that the distribution of F n converges to T u in the sense of the distance d^ . 



2.4 An extension of de Jong's theorem for degenerate U-statistics 

In the present and the subsequent section, we shall work within the following framework. We fix an 
integer d > 1, and let Y = {Yi : i > 1} be a sequence composed of i.i.d. random variables with values 
in M. d , whose common distribution has a density p(x) with respect to the Lebesgue measure on M. d 
(written dx). The sequence {N(n) : n > 1} of integer- valued random variables is independent of Y 
and such that, for every n, N(n) has a Poisson distribution with parameter n. It is well-known that, 
in this setting, the random point measure 

N(n) 

(2-8) V n := 

i=i 

(where 5 y represents the Dirac mass at y) is a Poisson measure on Z = M. d (equipped with the 
standard Borel a-field &(R) m ) with control measure /j, n (dx) = np(x)dx. We shall also use the 
shorthand notation //(dx) := fi\(dx) = p(x)dx. 

Our aim below is to provide a Gamma-type counterpart to a famous theorem by P. de Jong, proved 
in [3], involving sequences of degenerate [/-statistics of order 2. We stress that the results contained 
in [3] have later been extended to degenerate [/-statistics of a general order; see [U Albeit our 
method clearly applies to these general objects, we prefer here to focus on [/-statistics of order 2, in 
order to obtain neater statements and to emphasize the method over technical details. We start with 
some useful definitions. 

Definition 2.10 ([/-statistics). (i) Let k > 2, and let h : M. q — > R be a symmetric kernel such that 
h G Lg ym (/i fe ). The (symmetric) U -statistic of order k based on h and on the sample {Yi, . . . , Y m } 
(where m > k is some integer) is the random variable 

m 

(2.9) U m (h,Y)= h(Y il7 ...,Y ik ), 

»i,...,ifc=l 

where the symbol indicates that the sum is taken over all vectors . . . such that 
ij 7^ ig for every j ^ I. 

(ii) Fix k > 2 and let U m (h, Y) be a symmetric [/-statistic as in (|2.9[) . The Hoeffding rank of 
U m (h, Y) is the smallest integer 1 < q < k such that E[h(Y 1 , . . . , Y k )\Y u . . . , Y q _{\ = (a.s.-P) 
andE[/i(Yi,...,Y fc )|Yi,...,r,] / 0, where E[h(Y u Y k )\Y u Y ] := E[/t(Y 1; . . . , Y k )]. A U- 
statistic of order k with Hoeffding rank equal to k is said to be completely degenerate. In other 
words, a [/-statistic such as (|2.9p is completely degenerate if h is a non-zero kernel verifying 

h(x,yi,. . . ,y k _ 1 )p(x)dx = (^ fc_1 - a.e.). 

(iii) A collection of random variables {F n : n > 1} is said to be a sequence of geometric U -statistics 
of order k, if there exists a kernel h G Lg ym (/x fc ) such that 

F n = U N{n) (h,Y), n>l, 

where {N(n) : n > 1} is the independent Poisson sequence introduced above. 
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Before presenting the main result of this section, and in order to make the connection with our 
general framework more transparent, we shall recall an important finding from [30, Lemma 3.5 and 
Theorem 3.6], stating that Poissonized [/-statistics of order k live inside the sum of the first k + 1 
Wiener chaoses associated with the Poisson measure %. The proof heavily relies on results by Last 
and Penrose |12j . 

Lemma 2.11 (Reitzner and Schulte). Consider a kernel h E Lg ym (/x fc ) such that the corresponding 
Poissonized U -statistic U^f^(h,Y) is square-integrable. Then, h is necessarily in L 2 ym (^L k ), and 
C/jV(n) (h, Y) admits a chaotic representation of the type 

k 

U N{n) (h,Y) = E[U N{n) (h,Y)} + Y j n k - i I l (h t ) 

i=l 

where li indicates a multiple Wiener-Ito integral of order i with respect to the compensated Poisson 
measure fj n = rj n — fi n , defined according to \2. 8\) . and 

(2.10) ft* fa,... = (^\ J^ h{z l ,...,z i ,»)v k - i (d»), (z u ...,Zi) G Z\ 

where the bullet stands for a packet of k — i variables that are integrated with respect to [i k ~ l . In 
particular, h = hy- and the projection hi is in L 2 (/j, 1 ) for each 1 < i < k. 



The following statement corresponds to the main result proved by de Jong in [3], in the special case 
of symmetric U -statistics of order 2 (note that the assumption that the underlying kernels have finite 
moments of order four is only implicit in de Jong's work). Given positive sequences a n ,b n , n > 1, we 
write a n ~ b n whenever lim a n /b n = 1. 

n— >oo 

Theorem 2.12 (de Jong). Let {h n : n > 1} be a sequence of non-zero elements of L^ ym (fi 2 ) . Define 
Fn = U n (h n ,Y) and assume that F n is completely degenerate. Then, one has that a 2 (n) := Var(i ? „) ~ 
2n 2 ¥,[h n (Yi, Y%) 2 ], and the fourth moment condition 

n->oo a{np 

implies that, as n — > oo, the sequence F n := F n /a(n) converges in distribution to a standard Gaussian 
random variable. 

The following statement consists of two parts. Part (A) is a quantitative extension of Theorem 
I2.12l based on a direct study of the fourth moments of the Poissonized [/-statistic, whereas part (B) is a 
Gamma-type extension of de Jong's theorem which is directly based on the results discussed in Section 
12.11 Apart from [1], our findings should be compared with the seminal work by Jammalamadaka and 
Janson [8] , about the normal and Poisson approximation of U -statistics of order two. To our knowledge, 
the forthcoming Theorem 12.131 is the first quantitative extensions of the de Jong theorem, also dealing 
with the non- normal approximation of general degenerate [/-statistics. Moreover, we would like to 
emphasize that our proof of Part (A) is shorter and more transparent than the one presented in the 
original work [4] (one should note that, however, our methods only allow us to deal with symmetric U- 
statistics). Recall that the Wasserstein distance between the laws of two integrable random variables 
X, Y is given by 

d w (X,Y):= sup \M[h(X)]-M[h(Y)]\, 
/ieLip(i) 

where Lip(l) is the set of Lipschitz functions h : R — > M with a Lipschitz constant < 1. Recall that, 
in the framework of this section, Z = M. d . 
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Theorem 2.13 (Extended de Jong theorem). Let {h n : n > 1} be a sequence of non-zero elements 
o/L 4 ym (// 2 ) such that 

SUp 2 < °°- 

n (f z hi dfil) 

Put F n = U n (h n ,Y) and F' n = U ^ ( n )(h n ,^[) , and assume that these U -statistics are completely 
degenerate. Then, a(n) 2 := Var(i ? n ) ~ Var(i^) = 2n 2 K[h n (Yi, Y2) 2 ], and the following two points (A) 
and (B) hold. 

(A) // 

K\( F' ) 4 1 

(2.11) l \ 7. J -» as n^oo, 

cr(n) 4 

then both F n := F n /a(n) and F' n := F^/a(n) converge in distribution to a standard Gaussian 
random variable N. Moreover, there exists a universal finite constant K, independent ofn, such 
that, as n — > 00, 

(2.12) d w (K,N) <K x B n — > 0, 

(2.13) d w (F n ,N) <Kx (B n + n- 1 ^) — ► 0, 

with B n := a(n)~ 2 max | ( J 2 /i 4 d,u 2 ) 1/2 ; \\h n *{ h n \\; \\h n *\ h n \\\. 

(B) // J z hf l dfx 2 l —?• and i/iere exists v > suc/z f/iaf <r(n) 2 — >■ 2i/, and 

(2.14) E[(F^) 4 ] - 12E[(F^) 3 ] — ► 12^ 2 - 48^ as n -> 00, 

i/jen 6oi/i -F n and i 7 ^ converge in distribution to a random variable G{y), which has distribution 
T v . Moreover, there exists a universal constant K > such that, as n —> 00, 

(2.15) <k{K,T») < ciMK) + <%MK) * A 5 (F^)<KxC n — ► 0, 

(2.16) d 3 (^n,r,) < ^ x (C n + n" 1 / 4 ) — > 0, 

with C n := max||2||/i n || 2 — 2^|; ( ^nd/^n) 1 ^ 4 ^!^ *2 ^"I^^HI^n^i^n — we have 

used the notation introduced in \2. ip - JjQ)) . 

Remark 2.14. Our proof of Theorem 12.131 shows indeed that the quantity -B n (reps. C n ) in the 
statement converges to zero if and only if the asymptotic condition (|2.1ip (resp. (|2.14p ) is verified. 

2.5 Gamma convergence of geometric [/-statistics: characterization and bounds 

As anticipated, the aim of this section is to apply the main estimates of the present paper in order 
to characterize the class of geometric [/-statistics based on Y converging in distribution towards a 
Gamma random variable. Since our analysis is based on Theorem l2.13l our results will provide explicit 
estimates on the speed of convergence. We refer the reader to j6[ [32] for some classic references on 
the subject and to [11[ [30] for a discussion of several recent developments. We let the notation and 
assumptions of the previous section prevail and recall that a Gaussian measure G on (M. d , ^(R)® rf ) , 
with control fj,(dx) = p(x)dx, is a centred Gaussian family of the type 

G = {G(B) : B £ <%(R)® d , y,(E) < 00} 

such that, for every m > 1 and every B\, . . . , B m £ S3(M.)® d with n{Bi) < 00 (i = l,...,m), the 
vector (G(-Bi), . . . , G(B m )^ has an m-dimensional joint Gaussian distribution with covariance matrix 
E[G(Bi)G(Bj)] = n(Bi n Bj). 

The next statement combines findings from |11[ Section 7] (point (i)) with a classic characterization 
of elements in the second Wiener chaos of a Gaussian measure (point (ii); see [181 Section 2.7.4] for 
more details. 
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Proposition 2.15. Let k > 2 and let h G Ll ym (n k ) be a non-zero kernel such that the U -statistic 
F' n := U^f n )(h,Y) is square-integrable for every n, and has Hoeffding rank equal to 2. For n > 1, 
define also the standardized U -statistic F' n = n 1 ~ k F! n . 

(i) For every n, there exists a sequence of double integrals I 2 (fn) (each realized with respect to the 
compensated Poisson measure rj n — /j, n ) such that, as n — > oo, E[(J^ - h{fn)) 2 ] -> 0. Moreover, 
F' n converge in distribution to I^ih^), where L± indicates a double Wiener-Ito integral with 
respect to the Gaussian measure G, and h 2 is defined according to H2.10\) . The same convergence 
takes place for the de-P oissonized U -statistics F n = n 1 ~ k F n , where F n := U n (h, Y). 

(ii) The random variable (^-2) cannot be Gaussian. Moreover, assume that (^2) follows a 
T v - distribution. Then, necessarily, u G {1,2,...} and there exists an orthonormal system 

V 

{ei, . . . , e u } C L 2 (fi) such that h 2 = ej <g) ej and E[ej(Y"i)] = for i = 1, . . . , v. 

i=l 

V 

Remark 2.16. Let k = 2, and consider h.2 = Yl e^ ® e^, as in the statement of Proposition 12. 151 - (ii) ■ 

i=l 

Then, it is easily seen (by a direct computation) that 



U n (h 2 ,Y) 



V 

i=i 



n 



M{Y k )\ -}e t (Y k ) 



k=l 



The fact that the distributions of F' n and F n converge to T u is therefore a direct consequence of the 
usual multidimensional central limit theorem and of the law of large numbers. 

The next statement is a quantitative counterpart to Proposition 12. 15T - (ii) . containing in particular 
estimates involving Malliavin operators. Such estimates will be put into use in the Examples 12.221 - 
12.261 below, where the asymptotic behavior of a [/-statistic such as [//v( n )(^2, Y) is studied within the 
framework of hybrid convergence in random graphs, random flat and random simplex models. 

Theorem 2.17 (Bounds on Gamma convergence). Let the assumptions and notation of Pro- 
position \2A5\ prevail. Assume moreover that L 2 - ! (h2) has distribution T v for some v = 1,2,..., and 
also that {e\,...,e u } C L A (jx), where the orthonormal system {e\, . . . ,e u } is defined in Proposition 
\2.15\ -fii). Then, there exists a finite constant K, independent of n, such that 

d 3 (K,T u ) < cyAxiFi) + c 2 A A (F^) x As{F' n ) <Kx n~ l '\ 
d 3 (F n ,T„) < Kxn- l l\ 

where in the first inequality we used the notation defined in \2. l\) - [2~3\) 

Example 2.18. (i) Let 51,52 be two orthonormal elements of L 2 (/j.) such that 51,52 £ -^ 4 (m) arm 
f z gi(z)/j,(dz) = J z g2(z)n(dz). We stress that we do not require that 51,52 have disjoint 
supports. Then, the kernel 

hi \ l ( \^.( w \ -52(21) 51(22) -52(^2) 
h 2 (z 1: z 2 ) = - (51 - 52) <8> (51 - 52) (21, z 2 ) = y= x 



is such that the corresponding [/-statistics of order two F n := U n (h 2 , Y) and F' n := [/jy(n)(^2, Y) 
are completely degenerate, and both converge in distribution to Ti, with an upper bound of order 
n -1 / 4 on the rate of convergence. 

(ii) As an example of a pair (51, 52) verifying the requirements at Point (i), one can take 51 = \/2 1a 
and 52 = \/21_B) where {A,B} is a measurable partition of Z such that n{A) = /x(-B) = 1/2. 
Considering the case d = 1, p(-) = \ l(-i,i)(- ), 51(2) = V21(o,i)(^) and g 2 {z) = \f2 l(_i )0 )(^), 
one obtains a kernel hi with support in (—1, 1) 2 \{(0, 0)} and such that h 2 {z\, z 2 ) = 1 if z\z 2 > 0, 
and ^2(21, 22) = —1 if z\Z2 < 0. In this way, one recovers the non-central result discussed by 
Reitzner and Schulte in [301 en d of Section 5.1]. 
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2.6 Multivariate extensions and hybrid convergence 

We describe here three multivariate extensions of the results in Section 12.11 The first two results can 
be seen as partial analogues on the Poisson space of [191 Theorem 4.4] - for the multivariate Gamma 
convergence - and |19[ Theorem 4.5] - for the hybrid convergence - both concerning sequences of 
multiple integrals with respect to a Gaussian measure. Observe that the method used in [19] is based 
on new criteria for asymptotic independence of multiple integrals. These criteria are not available 
on the Poisson space. For this reason, our approach is different and will be in the spirit of the 
"interpolation method" used in [2J. Such an interpolation method will be also used to deduce our third 
result, concerning hybrid Poisson/Gamma convergence. As already pointed out, the general problem 
of characterizing independence is rather well-understood in a Gaussian framework (cf. l36| I37j). 
while the topic is still largely open in the context Poisson measures; see [28[ I31| . 

Remark 2.19. As before, we consider the framework of a sequence of Poisson measures {rj n : n > 1} 
(on some Polish space (Z,3f)), each having a finite with non-atomic control measure /x n . For the 
entire section, I q denotes the multiple Wiener-Ito integral, of order q, with respect to one of the 
compensated measures fj n = rj n — [A n (the concerned index n will always coincide with the index of the 
integrated function, for instance: I q (f n ) indicates the multiple integral of order q of f n with respect 
to ii n ). 

Let d > 1 be a fixed integer, let v\,...,v^ > and let (G?i, . . . , Gd) be a vector consisting of 
independent random variables such that G{ has the centred Gamma distribution T Ui . Further let 2 < 
Qi < 12 < ■ ■ ■ < Qd be even integers satisfying 2qi ^ qj for any i ^ j and let for % G {1, . . . , d}, {fn^ : n > 
1} be a sequence of kernels such that fjp G Lg ym (/Un ), satisfying in addition the technical conditions 
of Section [3}- (VIII). The next result deals with the announced multivariate Gamma convergence. We 
emphasize that we do not need further conditions on asymptotic covariances due to our assumption 
that all multiple integrals have different orders. 

Theorem 2.20 (Multivariate Gamma convergence). Let the above notation and conditions pre- 
vail. For any i 6 {1, . . . , d} assume that lim ^jlll/i, || 2 = 2z/j and that 

lim ** /« || = and lim ||/« Sjg/W - c q J^ || = with c qi = J 

\2)'\ qi /2) 

for all pairs (ri,£i) such that either = qi and li = 0, or r{ £ {1, . . . ,qi}, 1% 6 {1, ■ ■ ■ , min(rj,gj — 
1)} and Ti and li not equal to q-i/2 at the same time. Then, (iqxifn ),-■■, Iq d (fn)) converges in 
distribution to (G±, . . . , Gd) as n —> oo. 

We go one step further and turn to an extension of Theorem 12 . 201 where we consider convergence of 
a random vector of multiple integrals to a hybrid random vector whose components are independent 
and in part centered Gamma and in part standard Gaussian random variables. A similar setting with 
Poisson random variables instead of centred Gamma ones has recently been studied in [2]. However, 
we would like to emphasize that, in contrast to [2], here both distributions considered in the target 
vector are absolutely continuous with respect to the Lebesgue measure on the real line. 

To formulate our result on the Gamma/Gaussian hybrid convergence, let d\,d2 > 1 be fixed 
integers, let v%, . . . , > and let (G±, . . . , G^ , N^+i, • • • , Ndj^+d 2 ) be a vector consisting of inde- 
pendent random variables such that Gi has distribution T Ui for i 6 {1, . . . , d%\ and iVj has a standard 
Gaussian distribution for i G {d% + 1, . . . , d\ + d%\. Further let 2 < q±, q2, . . . , qdi+da be integers 
such that the following constraints are verified: (a) qi < ■ ■ ■ < qd 1 and, in general, qi ^ qj for every 
1 < i 7^ j < ^i + d,2, (b) qi, . . . , qd 1 are even integers, and (c) there is no pair (i, j) G {1, . . . , d\ + c?2} 2 
such that i € {1, . . . , d±} and 2qi = qj. For i G {1, . . . , d\ + c?2}, we let : n > 1} be a sequence of 

symmetric and square- integrable kernels such that f$ G Lg yra (^n), satisfying moreover the technical 
conditions stated in Section [3]- (VIII). 
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Theorem 2.21 (Gamma/Normal hybrid convergence). Let the above notation and conditions 
prevail. Assume that lim (?!||/n II 2 = 2^ whenever i G {l,...,di} and that lim qj'l 1 '""' 1 - 1 

n— >oo ' . — 

whenever i G {cq + 1, . . . , c?2}. Furthermore, for any i E {!,... , di} suppose that 



n— ¥oo 



(2.17) lim /« ||=o and Urn ||/« - c 9i /«|| = with c q , 

2 



/or aZZ pairs (ri,£i) such that either ri = qi and l{ = 0, or r» € {1, ... , a;}, 4 G {1, . . . , min(rj, a; — 1)} 
and and £i not equal to Oj/2 at the same time. For i G {d\ + 1, . . . , d\ + a^} suppose that 

(2-18) lim [|/»*S/W|| = 



for all ri G {1, . . . and^ G {1, . . . , min(r i; %-l)}. TTiera j ), . . ■Jg dl+d2 (fn 1 )) converges 



m distribution to [Gi, . . . , , A^+i, . . . , A^+da) as re — )■ 00. 

Example 2.22. We illustrate Theorem 1 2 . 2 1 1 wit h an example related to the theory of random graphs; 
the reader is referred to [25J for an introduction to this topic. Note that we will allow the underlying 
Poisson measure to depend on re; see Remark 12.191 Let d > 1, and define Y and r/„ as in Section f2.41 
Let {r n : n > 1} be a sequence of strictly positive numbers decreasing to zero. For every n, we define 
D n ■= (V n ,E n ) to be the random 'disk graph' obtained as follows: V n = {Yi : i = 1, . . . , N(n)} and 
two vertices Yi,Yj G V n are connected by an edge if and only if their Euclidean distance is strictly 
positive and less than r n (in particular, D n has no loops). Now let A be a feasible connected graph 
(in the sense of [21 [TT| I25j ) with q vertices, where g / 2,4. For every n, we define L n to be the 
random variable equal to the number of induced subgraphs of D n that are isomorphic to A, that is, 
L n is equal to the number of subsets of the type Yr q \ = {Y^, . . . , Yi q } C V n such that the restriction 
of D n to is isomorphic to A. We also set L n = (L n — E[L n ])/Var(L n ) 1//2 . Now assume that 
nr^ — > and n q (rfy q ~ l -4 00. According to the discussion contained e.g. in [251 Chapter 3] or 
Section 3], one has that the following four facts are in order: (i) E[L n ] ~ KQn q (r^) q ~ 1 (for some finite 
constant Kq > 0), (ii) Var(L n ) Kin q {r a \) q ~ l (for some finite constant K\ > 0), (iii) there exists 
a sequence of multiple integrals of order q with respect to fj n = r\ n — u n , say I q (f n ), such that, as 
n — >• 00, E[(L n — Iq{fn)) 2 ] — > 0, and (iv) the kernels {f n } verify the asymptotic relation (|2.18p (where, 
for every n, the contractions and norms have to be considered with respect to the measure fj, n ), so 
that L n converges in distribution to a standard Gaussian random variable as n — > 00. Now consider a 
sequence {F^} of degenerate [/-statistics of order 2 as in Theorem 12.171 (for instance, those appearing 
in Example l2.18p . Since each F' n is a double Wiener-Ito integral (with respect to r] n ) verifying condition 
(|2.17p and q ^ 4, we can directly apply Theorem l2.21l in the case d\ = o\ = 1, q\ = 2 and qi = q, and 
conclude that, as n — > 00, the pair (F^,L n ) converges in distribution to a vector (G,N) composed 
of independent random variables such that G has distribution T u and N follows a standard Gaussian 
random variable. 

We finally show how one can use the results of the present paper to deal with a hybrid Pois- 
son/Gamma convergence (we just consider two-dimensional vectors in order to simplify the discussion, 
but there is no additional difficulty in considering vectors of higher dimensions) . Let v, A > and let 
(G, P) be a vector consisting of independent random variables such that G has distribution T u and 
P has a Poisson distribution with mean A. We fix an even integer q > 2, and consider a sequence 
{fn '■ n > 1} of kernels with f n G Lg ym (^n) and such that the technical conditions stated in Section 
[3}-(VIII) are satisfied. We also consider a sequence {H n : n > 1} of random variables such that: (a) 
each H n is a functional of the Poisson measure rj, which is in the domain of the Malliavin derivative 
D and takes values in Z + = {0, 1, 2, . . .}, (b) the numerical sequence 

(2.19) n^E J^(D z H n ) 2 fi n (dz), n > 1, 

is bounded. 
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Theorem 2.23 (Gamma/Poisson hybrid convergence). Assume that lim g!||/n|| 2 = 2f and 

n— >-oo 

(2.20) lim || / n ^ /nil = and lim \\f n ^J/lfn ~ CqfnW = with c q - 



/or all pairs (r, £) suc/i that either r = q and £ = 0, or r G {1, . . . , g} ; ^ G {1, . . . , min(r, q — 1)} and r 
and £ not egua/ tog/2 at i/ie same time. We also assume that, as n —> oo, 



(2.21 



E[# n ]^A, E |A — (DH n , —DL~ X H^)\ ->■ 0, 

and E^ \D z H n (D z H n - l)D z L~ l H n \ fi n (dz) -> 0. 

TTien (l q (f n ), H r ^) converges in distribution to (G, P) as n 



oo. 



Example 2.24. We consider the same framework and notation as in Example 12.221 Here, we take 
g > 2 to be a general integer (which can be possibly equal to 2 or 4), whereas A is a feasible connected 
graph of order q. We stress that, for every n, the random variable L n is a functional of the Poisson 
measure rj n on whose control measure is given by fx n {dx) = np(x)dx. We put r n = n~ dq '( q ~ 1 ', 
in such a way that n q (r d ) q ~ 1 = 1. According e.g. to the analysis contained in [251 Chapter 3] or [2j 
Section 2.4], one has that the following two facts are in order: (i) there exists a constant A > such 
that E[L n ] « Var(L n ) m Xn q (r d ) q ^ 1 = A, and (ii) the sequence H n = L n satisfies (|2.19p . as well as the 
asymptotic relations (|2.2ip (here, for every n, the Malliavin operators are defined with respect to the 
ran dom measure r\ n and the inner products and integrals are obtained by integrating with respect to 
M = Mn)> so that L n converges in distribution to a Poisson random variable with mean A. Considering 
a sequence {F^} of degenerate [/-statistics of order 2 as in Theorem 12.171 (see e.g. Example 12. 18f) . one 
has that each F' n is a double integral verifying condition (|2.20p . We can therefore apply Theorem l2.23l 
and infer that, as n — > oo, the pair (F^,L n ) converges in distribution to a vector (G,P) composed 
of independent random variables such that G is distributed according to T u and P has a Poisson 
distribution with mean A. 

Example 2.25. Let us consider a Poisson measure rj n of fc-dimensional flats in M. d with 2k < d (where 
the flats are suitably parameterized to fit into our framework) . We assume that the distribution of rj n 
is invariant under rigid motions for each n, and that rj n has intensity n > 1. Let us fix a closed convex 
set W C M. d with volume one and define the distance distw{E, F) of two fc-flat E, F as the minimum 
over the Euclidean distances of xe G E n W and xj? G Fn W . By M n we denote the number of 
pairs (E,F) of distinct flats of r\ n such that distp^(i£, F) < r n , where r n = n - 2 /( d ~ 2k ) . According to 
Theorem 2.1 in [31] we know that (i) there exists a constant < A < oo (depending on d, k and W) 
such that E[M„] ~ YM n ~ A, (ii) the asymptotic relations (|2.2ip are satisfied, and (iii) M n fulfills the 
technical condition (|2,19p . Thus, M n converges in distribution to a Poisson random variable with mean 
A. Let now {F^} be a sequence of degenerate [/-statistics of order 2 as in Theorem 12. 171 see Example 
12.181 Then, as in the previous example, each F' n is a double integral such that condition (|2.20p is 
verified. Thus, Theorem 12.231 can be applied to show that the random vector (F^,M n ) converges in 
distribution to a random vector (G, P) with independent components such that G has distribution Y v 
and P has a Poisson distribution with mean A. 

Example 2.26. Let W C R d be a closed convex set with volume one and let Y be a sequence of 
i.i.d. points in W, which are uniformly distributed and whose random number N(n) follows a Poisson 
distribution with parameter n G N. Any d + 1 distinct points of Y form a non-degenerate random 
simplex in W . Define r n := n~^ d+1 ^ and let V n be the total number of such simplices whose volume does 
not exceed r n . Then (i) there exists < A < oo (depending on d and W) such that M[V n ] ~ YV n ~ A, 
(ii) the asymptotic relations (|2.21|> are satisfied, and (iii) V n fulfills the technical condition (|2.19p so 
that the law of V n converges, as n — > oo to a Poisson distribution with mean A. This can be seen from 
Theorem 2.5 in [34] . Define the sequence {F^} of degenerate [/-statistics as in Example 12. 181 or. more 
generally, as in Theorem 12.171 Then, following the same line of reasoning as as above, Theorem 12.231 
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can be applied to show that the random vector (F^, V n ) converges in distribution to a random vector 
(G, P) with independent components such that G has distribution T\ and P has a Poisson distribution 
with mean A. 

3 Background material 

In this section we collect definitions and results that are needed in the statements and proofs of our 
results. For more details, we refer to the monographs [231 EZ] or to the papers [T2"| [20]. 

(I) Poisson measures. We shall denote by i] a Poisson measure with non-atomic and u-finite control 
measure fx on some Polish space Z (which is endowed with the Borel u-field 3f). Recall that rj is a 
collection {17(B) : B G 3?$} of random variables indexed by the members of 2?o = {B G 3f : fJ-(B) < 00} 
such that: (a) i](B) follows a Poisson distribution with mean 71(B) for all B G 3fo, and (b) whenever 
A, B G 3?q are disjoint, i](A) and t](B) are independent random variables. By P^ we will denote the 
distribution of 77 (on the space of cr-finite counting measures on Z) . 

(II) L 2 -spaces. For q > 1 we denote by L 2 (fi q ) the L 2 -space L 2 (Z q ,^® q ,[x q ) and by L 2 yia (fi q ) 
the subspace of L 2 ([i q ) consisting of functions that are fi q -&.e. invariant under permutations of its 
arguments, so called symmetric functions. Suppressing the dependency on q, the scalar product and 
the norm in L 2 (/j, q ) (and L 2 ym (fi q )) are denoted by (•, •) and || • ||, respectively. In addition, we let 
L 2 (P V ) be the space of square-integrable functionals of 77. To avoid confusion we will use capitals to 
indicate elements of L 2 (P V ) and lower cases for elements of L 2 (/j <? ) or L 2 ym (/j, q ). We finally introduce 
the space L 2 (P, L 2 (/x)) = L 2 (J7 x Z^T ® Peg) 71) as the space of jointly measurable mappings 
u : x Z — > R such that E f z u(z) 2 fi(dz) < 00 (recall that (0,,F, P) is the underlying probability 
space). 

(III) Multiple stochastic integrals. For every integer q > 1 and every deterministic function 
/ G L 2 ym (/j, q ) let us indicate by I q (f) the multiple Wiener-Ito stochastic integral of order q of / with 
respect to the compensated Poisson measure 77 — 71. For general / G L 2 ([i q ) we put I q (f) := Iq(f), 
where f(x±, . . . , x q ) = (<?!) Ylm f( x n(i)i ■ ■ ■ > x n(q)) is the canonical symmetrization of / and the sum 
in its definition runs over all q\ permutations tt of {1, . . . , q}. The multiple stochastic integrals satisfy 
the following properties: 

Ef/J/i)] = (i = 1,2), and E[l qi (f 1 )I q , 2 (f 2 )] = q\\(h,h) l(<?i = 92) 

for any q u q 2 G {1,2, . . .}, h G L 2 ym (fi q ^) and f 2 G L s 2 ym (/j^). 

(IV) Chaotic representation property. The q-th. Wiener chaos W q associated with the Poisson 
measure 17 is the Hilbert space W q = {I q (f) : / G L 2 yra (fi q )}. In addition, we put Wo ■= M. It is 

a crucial property of 17 that L 2 (P V ) can be written as a direct sum of Wiener chaoses, i.e. L 2 (i-* r? ) = 

00 

® Wq. As a consequence, every F G L 2 (P V ) admits a (unique) chaotic decomposition in the sense 

9=0 

that 

00 

(3-1) F = E[F] + y £lM 

q=l 

with suitable functions f q G L 2 ym (/j <? ) and where the series converges in L 2 (P V ). 
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(V) The Malliavin operators D, L 1 and 5. The domain domD of the derivative operator D 

oo 

is the set of all F G L 2 (P V ) admitting a chaos decomposition (|3.ip such that ^ 9 < °°- For 

such F the random function Z 3 z ^ D Z F is defined by 

oo 

(3.2) £> z F = 5>/,-i (/,(*, •)). 

9=1 

where • ) is the function f q with one of its argument fixed to be z. Notice that DF G L 2 (P, L 2 (/j,)). 
The derivative operator can be also characterized as an "add-one cost operator", as follows; see [12|, I20J 
for proofs of this fact. For F G L 2 (P 7? ) and let F z (rj) be the random variable F(n + Then, 

for F G domL> and /u-almost every z £ Z, we have the identity .D z -F = F z — F, a.s.-P. Throughout the 
text, we also implicitly use the following converse statement (the proof is an elementary consequence 
of the main findings of [12] , and is included for the sake of completeness) . 

Lemma 3.1. Let F G L 2 (P V ) be such that E J Z (F Z - F) 2 n(dz) < oo. Then, F G domD. 

Proof. For every z G Z, define the 'trajectoriaF difference operator D' z F(n) = F z {rf) — F{rf). According 
to |X2^ Theorem 1.3], the square-integrable random variable F admits a chaotic decomposition of the 
type ()3.ip . with 

f q (z 1 ,...,z q ) = ±E[D' Zi ...D' Zq F], q = l,2,... 

(in particular, the deterministic function on the right-hand side of the previous equation is a well- 
defined element of L 2 ym (^ 9 ) for every F G L 2 (P V ) and every q > 1). In view of the assumptions, there 
exists a measurable set Z' such that \i{Z\Z') = and E[(D' Z F) 2 } = E[(F Z -F) 2 ] < oo for every z G Z'. 
It follows that the statement is proved once we show that, for every z G Z', the chaotic decomposition 
of D' Z F coincides with the right-hand side of f|3.2[) . Again by virtue of |12t Theorem 1.3], one has that 
the qth integrand in the chaotic decomposition of D' Z F is given by the mapping 

(*!, . . . , z q ) ^ ^ E[D' Z1 . . . D' Zq D' z F] = (q + l)f q+1 (z, Zl ,...,z q ), 

which yields the desired conclusion. □ 
For any F G L 2 (P V ) with chaotic decomposition (|3.ip satisfying E[F] = we put 

oo 

L- 1 F = -Y d Q- X I q {f q )- 

The operator L^ 1 is the so-called pseudo-inverse of the Ornstein- Uhlenbeck generator. Finally, we 
observe that, due to the chaotic representation property of f), every random function u G L 2 (P, L 2 (fj,)) 
admits a (unique) representation of the type 

oo 

(3.3) uz = 52l q (f q {z,')), zeZ, 

q=0 

where, for every z, the kernel f q (z, • ) is an element of L 2 ym (/i 9 ). The domain of the divergence 
operator, denoted by dom<5, is defined as the collections of those u G L 2 (P, L 2 (/i)) such that the 
chaotic expansion (|3.3p verifies the condition 

oo 

^( 9 + l)!||/ 9 || 2 <oo. 

g=0 
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If u G dom 5, then the random variable 8(u) is defined as 

oo 

<y(«) = X)i g+ i(/ g ) > 

q=0 

where f q stands for the canonical symmetrization of f q (as a function in q + 1 variables) . The following 
classic result, proved e.g. in |20j, yields a characterization of 8 as the adjoint of the derivative D. 

Lemma 3.2 (Integration by parts formula). For every G G domD and every u G dom<5, one has 
that 

(3.4) E[G<y(u)] = E[{DG,u)\, 
where, more explicitly, 

(DG,u) = J D Z G x u(z) fi(dz). 

(VI) Contractions. Let / G Lg ym (fj, q ) for some integer (? > 1 and r G {0, . . . ,p}, I G {1, ... ,r}. 
The contraction kernel on z 2q ^ r ~^ acts on the tensor product /<8>/ first by identifying r variables 
and then integrating out £ among them. More formally, 

/ *r /(7l> • • • ,Kr-lih, ■ ■ ■ ,tq-r,Sl, • • • , Sq- r ) = / f(zi, . . . , Z£, Jl, . . . ,7 r _£,tl, . . . , t q - r ) 

x f(zi, ■ ■ ■ , Hi 7i> • • • > SV-f , si, ... , s g _ r ) //(d(zi, . . . , ^)) . 

In addition, we put 

/ *° /(7i> ■ ■ ■ > 7n *1, • ■ ■ 5 si, ... , s g _ r ) = /(71, . . . ,7 r , ti, . . . , t g _ r )/(7i, . . . , 7 r , si, . . . , s q - r ). 

Besides the contraction / -kf, /, we will also deal with its canonical symmetrization /* r /, which is 
defined as 

(/ *r/)( X l> • • • > X 2q-r-i) = - r - £)\ E^ *^ /)( X 7r(l)' ■ ' ' > X 7r(2 9 -r-£) ) i 

where the sum runs over all (2q — r — t)\ permutations of the set {1, . . . , 2q — r — £}. 

(VII) Product formula. Let qi, q% > 1 be integers, fi G Lg ym (/u' ?1 ) and /2 G L^ ym (fL Q2 ) be as in the 
previous paragraph. In terms of the contractions of f\ and /2 one can express the product of I qi (fi) 
and Iq 2 {n) as follows: 

min(gi,g 2 ) / \ / \ r / \ 

(3.5) Iq,{h)Iq.{h)= £ r \ 9 r j^)EQW-M(/l^/ 2 ); 

see [231 Proposition 6.5.1]. In the particular case q\ = qi =: q and f\ = /2, we may define Gq/ := 
g!||/|| 2 and 

q r 



Glf:=J2J2l(2q-r-£ = p)rl( q ) Q/*f/ 

, — n o— n \ / \ / 



r=0 e=o 

for p G {1, ... j 2q}, which allows us to re- write (|3.5p in the more compact form 

2q 

(3-6) JJ(/) = £ J p(G?/)- 

p=0 
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(VIII) Technical assumptions. Whenever we deal with a multiple stochastic integral or a sequence 
Iqifn) of such integrals with /„ G L 2 ym ([in) we will (implicitly) assume that the following technical 
conditions are satisfied: 

i) for any r G {1, . . . , q}, the contraction f n -k q r ~ T f n is an element of L 2 (fi r n ); 

ii) for any r G {l,...,q}, I G {l,...,r} and {z\, ■ ■ ■ , z^q-r-e) G Z 2q ~ r ~ l we have that (|/ n | ~k~ 
\f n \)(zi, . . . , Z2q- r -t) is well defined and finite; 

iii) for any k G {0, . . . , 2(q — 1)} and any r and ^ satisfying k = 2(q — 1) — r — i we have that 



(/nO*, /n(z> ' )) 2 ^n(dz) < OO. 

We remark that (iii) is automatically satisfied if the control measure /i is finite (which is the case 
in our Examples 12.221 and I2.24ll2.26p . Intuitively, conditions (i)-(iii) ensure that every manipulation 
involving contraction kernels performed below is justified and is in fact valid. For the detailed role of 
these conditions and their implications we refer to [TO] or [22] . 

4 Proofs of the results 
4.1 Proof of Theorem I2T1 

Before entering the details of the proof of Theorem 12.11 we recall some facts related to Stein's method 
for the Gamma distribution established by Luk in [T3] ; see also Pickett [26J for refinements in the case 
of an integer- valued parameter v. We start by considering the second- order Stein equation 

(4.1) h(x-u)-E[h(G*(u))) =2xg"(x)-{x-v)g'{x), x > 0, 

where G*{v) = G(y) + v, with G{y) distributed according to T v , and h G H 3 . It is shown in [T4"l 
Theorem 1] that flU) admits a solution V h such that HV^Hdo < fll^ (i) l|oo for j = 1,2,3. Note that 
the assumption h G Ti 3 automatically yields that h has sub-exponential growth, so that [14] Theorem 
1] can directly be applied. 

Now, we turn to the first-order Stein operator T for T u , which acts on differentiable functions 
/ : R -> R. It is given by 

Tf(x) = 2(x + v)+f'{x) - xf(x), xGR. 

The associated first-order Stein equation is 

h(x) -E[h(G(u))]=Tf(x), xGR, 

where h G H 3 . For such an h, a solution Uh of the Stein equation - in what follows, sometimes called 
a Stein solution - is provided by 



(4.2) U h (x) 



-l{h{x) -E[h(G(y)))] : x < -v 
V^(x + v) : x > — v. 



Recall that the probability density of G{y) is given by 

9u{x) = ^^(x + v yftr-i e -W i {x> _ v} . 

Since for our choice of the test function h the mapping x h-> Uh(x) = V^(x + v) is bounded on (v, oo), 
we can use [351 Lemma 4] to deduce that, necessarily, 

U h (x) = \ / (h(y)-E[h(G(u))])g u (y)dy, x > -u, 

z[x + i>)g v \x) j _ v 
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yielding that x h-> Uh(x) is continuous on R, as deduced from a simple application of de l'Hopital's 
rule at x = —v. Also, £7^ is twice differentiable on R \ {^} and satisfies the estimates ||E/ft||oo — 
max(2,2/V) =: c , ||^||oo < max(l, 1/v + 1/u 2 ) = a and ||^'||oo < max(2/3, 2/(3u) - 3/V + 4/V 3 ) = 
C2 (here, c\ and C2 are the constants from Theorem 12. ip . We stress that, albeit Uh(x) is continuous on 
R, such a function is in general not differentiable at x = —v (it is however right- and left-differentiable 
at such a point). We remark that the quantities 2, 1 and 2/3 appearing in the constants co,c\ and 
C2 come from smoothness estimates for (|4.2p on the interval (— v, oo), whereas the presence of the 
constants 2/V, \jv + and 2/(3^) — 3/V 2 + 4/V 3 is explained by elementary estimates of (|4.2|) on 
the interval (— oo,— v\. 

Let now J 72 be the space of continuous functions / on R, which are twice differentiable on R\ {y\ 
and satisfy 

ll/lloo < c , H/'lloo < ci and [|/"||oo < P2- 
In the light of the previous discussion, we conclude that 

(4.3) d 3 (F,T u ) < sup \E[2(F + v)+f'(F) - Ff(F)] | , 

where F and v are as in the statement of Theorem 12.11 and where here and below f'(—v) stands for the 
left-sided derivative of / at —v, i.e. f'{—v) = hni f'{x). We also refer the reader to Lemma 1.3 in 

x /*—v 

[16J and the references cited therein. The estimate (j4.3[) is the starting point of the proof of Theorem 

EU 

Proof of Theorem \2.1\ We have to show that the right-hand side of (|4.3p is bounded from above by 
the right-hand side of (|2.ip . This is done by borrowing some ideas from [33J. To start with, consider 
/ € J- 2 , write F(rj) instead of F to emphasize the dependency of F on r\ and fix z G Z. Because 
of the non-differentiability of the Stein-solution at — u, we will have to distinguish the three cases a) 
F(rj) < -v, F(rj + 5 Z ) < -v or F(ij) > -u, F(r) + S z ) > -v, b) F(rj) < -v < F(rj + 8 Z ) and c) 
F(jj + 8 Z ) < — v < F{rj). For a) we use a Taylor expansion to see that 

D z f(F( V ))=f(F(r) + 6 z ))-f{F(r,)) 

= f'{F(r,)) (F( V + S z ) - F(n)) + R{F( V + 5 Z ) - F( V )) 
= f(F{r 1 ))D z F{r 1 )+R a {D z F(r ] )), 

where the reminder R a is such that |i? (x)| < ^H/"!^^ 2 = \c2X 2 \ recall that / is differentiable on 
R\{— f}, as well as right- and left-differentiable at x = —v. For case b) we also use a Taylor expansion 
to see that 

D z f(F(rj)) = f(F( V + S z )) - f(F(rj)) = f{F( V + 5 Z )) - f{-u) + f{-u) - f(F( V )) 
= f(—v+) (F(r) + 5 Z ) + u) + i/"(F^)) {F( v + 5 Z ) + vf 

+ f'(F(r,)) {-v- F^)) + \f{F^)) mf 
= f(F(r,))D z F( V ) - f'{F( V )) (F(rj + 5 Z ) + u) + f'(-v+)(F( V + S z ) + v) 

+ \f"{F{v)) {F(V + 8 z ) + uf + \f'{F{rf)) (F( V ) + vf 
= : f(F(rj))D z F{rj) + R b {F( V ), z,u) 

with some F G ( — v,F(r] + 5 Z )),F G (F(r/),— ^) and where f'(—v+) stands for the right-sided 
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derivative of / at —v. Similarly, in case c) we find that 

D z f{F(rj)) = f(F( V + 6,)) - f(F(r,)) = f{F( v + 5 Z )) - f{-u) + f{-u) - f(F( V )) 
= f'(-u-)(F(r, + 5 Z ) + u) + \f n {F{fj)) {F(rj + 5 Z ) + vf 

+ f'(F(r,)) ( - v - F( V )) + \f" {Fir,)) {-v- F(r?)) 2 
= f'{F{r,))D z F{r,) - f'[F(r,)) [F( V + 5 Z ) + u) + f '(-*-) [F(r, + 5 Z ) + u) 

+ \f"{F{ri)) (F(V + S z ) + uf + \f'{F{rjj) (F( V ) + vf 
=: f'(F(r,))D z F(r)) + R c (F(i]),z, u) 

again with some F € (F(rj + 5 Z ), —v) ,F€. ( — u, F(rj)) and where /'(— v— ) stands for the left-sided 
derivative of / at v. Summarizing, we conclude that 

(4.4) D z f(F( V )) = f(F( V ))D z F( V ) + R{F( V ),z, u) 

(recall that f'{—v) = f'(—v—) by convention), where the global reminder term R(F(t]), z,iA is given 
by 

R(F(r]),Z, v) = R a (F(l])) {1{ F ( V ),F( V +S z )>~u} + 1 {F{ V ),F{ V +8 z )<-u}) 

+ R b (F(r]),z,u) 1{f( v )<-u<f{ v +8 z )} + Rc{F{rj),z,is) 1 {f{v+Sz )<- u <f( v )}- 

We have seen that R a has the property that |i? a (x)| < ^||/"||oo x 2 = \c 2 x 2 . For Rf, and R c we notice 
that in these cases \F(rj + 5 Z ) + v\ < \D z F(r))\ and \F(i]) u\ < \D z F(rj) |, which together with the 
properties of / £ T 2 leads to the bound 

\R{F{rj),z,v) \ < ljC 2 \D z F(ri)\ 2 (l {F{v)iF{v+Sz)> _ u} + 1{f(t,),f( v +s z )<-u}) 
+ (2ci\D z F(r])\ +c 2 \D z F(r])\ 2 ) 1{f( v )<-,s<f( v +s z )} 
+ (2a\D z F(ri)\ +c 2 \D z F(r])\ ) 1{f(t,+S z )<~u<f( v )} 

< c 2 \D z F{rj) | 2 + 2ci|^-F(?7)|(l{ F ( J? )<_ iy<F ( ?7+< 5 z )} + ^-{F(r 1 +5 z )<-u<F( v )}) 

= c 2 1 D Z F( V ) | 2 + 2 Cl (£>,1 { f(„)>-„}) {D z F(rf)) . 

Using now the integration by parts formula from Malliavin calculus, fj3.4[) in Lemma 13.21 and simpli- 
fying the resulting expression we find 

E[Ff(F)) = E[LL- 1 Ff(F)] = E[ — 5(DL^ 1 F)f(F)] = E[(Df(F), —DL~ X F)] , 
which in view of (|4.4|) leads to 

E[(Df(F),-DL- 1 F)] =E[f'(F)(DF,-DL- 1 F)] + E[(R(F, z, v), —DL~ x F)\ . 
Consequently, because of the above estimate on \R[F(rj), z,v)\, 
\E[2(F + v) + f'(F)-Ff(F)]\ 

< \E[f'(F)(2(F + u)+ - (DF, -DL- l F))] | + \E[(R(F, z, v), -DL- l F)] \ 

< c{K\2(F + u) + - (DF,-DL- 1 F) \ + c 2 j E\\D z F\ 2 \DL~ l F{\ fi(dz) 

+ 2 C1 J E [{D z l {F> _ v} ) {D z F)\D z Lr x F\] n(dz). 

This shows the first inequality (|2,ip in Theorem 12. 11 The second estimate ()2.2p follows from (|2.ip and 
the assumption that E[{DF, -DL~ l F)\F] > 0. This proves Theorem EH □ 
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4.2 Proof of Proposition 12.31 

We start by observing that the function x i— )■ 3>(x) := = sign(x) x 2 , x G R, is such that, for every 
a,b£R, $(6) = ${a) + 2\a\{b-a) + R{a,b), where |/2(o,6)| < (b-a) 2 . It follows that 

(4.5) ($(6) - $(a)) 2 < 8a 2 (6 - a) 2 + 2(6 - a) 4 . 

Since /i is finite, 

E / {D z l {F> _ v} f fi(dz) < E / (l {F+1)zF> _ i , } - l {F> _ !/} ) 2 ^(dz) < < oo, 

Jz Jz 

which implies that \{ F> _ U } G domD; see Lemma 13.11 and compare with Remark 12.21 (i). Moreover, 
our assumptions imply that DF\DF\ = &(DF) E dom<5. We can now apply the integration by parts 
formula (|3.4p . together with the relation L~ 1 F = —q~ l F, to deduce that 

qxA 3 (F) = eJ (D 2 1{f>-v})$(D z F) ii(dz) 
= E[1 {F> _ V} 5($(DF))) 
< [K[5{^(DF)) 2 }] 1/2 . 

Again in view of our assumptions, the Skorohod isometry implied by [27, Proposition 6.5.4] is verified, 
and we deduce that 

E[6{${DF)) 2 ] < E [ &(D Z F) 2 fi(dz) + E j I ' [D Z2 <S>(D Zl F)} 2 / u(dz 1 )/i(dz 2 ) 
= E f (D z Ffn(dz)+E f j [D Z .MD Z1 F)} 2 ^dz 1 ) f i(dz 2 ). 

Since D Z2 $(D Zl F) = $(D Zl F + D Z2 D Zl F) - ®{D Z1 F), we can now apply (g3]) with a = D Z1 F and 
6 = D Zl F + D Z2 D Zl F to infer the upper bound 

[Z^OD^F)] 2 < 8(D Zl F) 2 (D Z2 D Zl F) 2 + 2(D Z2 D Z1 F)\ 

and the conclusion follows immediately. □ 

4.3 Proof of Theorem E32 

Let F n = /,(/„) be as in the statement of Theorem [221 Then (DF n , —DL~ l F n ) = ±\\DI q (f Tl )\\ 2 and 
E[(DF n , -DL- 1 F n )\F n ] > 0. Thus, we need to prove that for such F n the right-hand side of QggD 
converges to zero as n — > oo. We do this by showing that the three terms Ai{F n ), A^{F n ) and A±{F n ) 
(see (|2.3p ) all converge to zero as n — > oo; the computations performed below will also implicitly 
provide the upper bound (|2.6p . It is important to note that our analysis of the terms Ai(F n ) and 
^4-4 {Fn) does not make use of the fact that fj, n {2) < oo. It is convenient to start with the reminder 
term A 4 {F n ). 

Lemma 4.1. Under the conditions of Theorem \2.(A it holds that A±{I q {f n )) — > 0, as n — > oo. 
Proof. First observe that in our case 



\ 



E[\D z I q {f n )\*] Mn (dz). 



z 



We can now use [22j formulae (4.17) and (4.18)] to deduce that 

q r — 1 

MWn)) < ^ ^ 1(1 < r + ^ < 2g - 1) 
(4.6) r=1 ^=o 



((r+/-i)o 1/, h-«-i)'( I iTl < )'( , ;l;')i/.*j/. 
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Since this estimate does not involve the middle contraction fn^M fni the- conclusion follows immedi- 
ately. □ 

Now we study the convergence of the sequence A'i(F n ). 

Lemma 4.2. Under the conditions of Theorem \2.6\ we have A' 1 (I q (f n )) — > 0, as n — > oo. 

Proof. One must prove that E[||D/ q (/ n )|| 2 — 2ql q {f n ) — 2qu\ 2 — > 0. Expanding the square and using 
the fact that E[I q (f n )] = we have to show that 

(4.7) E[\\DI q (f n )f] - AqE[I q (f n )\\DI q (f n )\\ 2 } + 4<z 2 E[I 2 (/ n )] - AqvE[\\DI q (f n )\\ 2 } + 4gV 

as n — > oo. Firstly, E[I 2 (f n )\ = ql\\f n \\ 2 -> 2v. The definition of DI q (f n ) and formula (J3SD imply that 

2(9-1) 



(4.8) 



\DI q (f n )\\ 2 =qq\\\f n f + q 2 E I I^Gf 1 f n (z, • )) Mn(d* 

p=l ^ 2 



so that E[||L'/q(/ n )|| 2 ] = <7 glU/nll 2 , which asymptotically behaves like 2gi/. Using integration by parts 
([33D together with the relation DF 2 = 2FDF + (DF) 2 applied to F = I q (f), we infer that 

E[I q (f n )\\DI q (f n )\\ 2 } = q -E[ll(f n )\ - l -El D z I 3 q (f n )vn(dz). 

Now, in view of the estimate (|4,6p . the second summand on the right-hand side of the previous 
equation converges to zero as n — > oo, and consequently E[l g (/ n )||L>/ g (/ n )|| 2 ] behaves asymptotically 
as |E[l| (/„)]. Using (|3.5p and the orthogonality of chaoses we obtain 

2 P 



P =o w ^= w 

^ ! C) 2 (A) 9!W,/n) ' 



p=q/2 

so that E[/ g (/„)||Z)/ g (/ ra )|| 2 ] has the same limit as 



p=q/2 

Moreover, one can show that 

2( g -l) 

(4-9) E[\\DI q (f n )f] = q 2 (ql\\f n \\ 2 ) 2 + q 4 E P«ll^/»ll 2 . 

P =l 

where Gpf n with p 6 {1, . . . , 2{q — 1)} is defined by 

q min(t,g-l) , _1\ 2 //— 1\ 

= £ E l(2(Z-*-fl=p)(t-l)! (?_ ) (.J/n^/n. 

t=i s=i V / V / 

Indeed, use (|4.8|) . the orthogonality of the random variables 

Ipi (G^Uz, • )) Hn(dz) and / J P2 {G^Uz, • )) Mn (dz) 
2 Jz 
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for 1 < pi 7^ p2 < 2(q — 1) as well as the stochastic Fubini theorem [23\ Theorem 5.13.1] (which is 
valid thanks to our technical assumptions made in Section [3|) to conclude that the identity (|4.9p is 
verified; see also the proof of Theorem 4.2 in [22]. We now exploit the assumption that ||/ n *r fn\\ — > 
with r and £ as in the statement of Theorem 12.61 It implies that 

(4.10) (/n*f P /n,/n) ^0 and (fn* S tfn,fn*(>fn)^0 

as n — > oo for all p G {q/2 + I, . . . ,q} and t, t' G {l,...,q}, s G {1, . . . , min(t, q — 1)}, s' G 
{1, . . . , min(t', q — 1)} and t, s, t 1 ', s' not equal to g/2 at the same time. Indeed, 

K/n*r P /-/")l ^ ll/n^" P /n|| ||/n|| < H/n^/nH ||/»|| 

for p G {g/2 + 1, . . . , q} and similarly 

\(fn* S tfn,fn* S t'fn)\ < [|/n*t/n|| ll/n *t<7n[| < [|/n*?/n|| [|/n*£ /n|| "> 0, 

where t,s,t',s' are as above. Plugging the expressions for E[||D7 9 (/ n )||], E[/ g (/ ri )||D/g(/ n )|| 2 ] and 
E[||D/g(/ n )|| 2 ] into ([4.7p and using the first statement in (|4.10p we see immediately that (|4.7p has the 
same limit as 

f x2 2(g-l) 
(4.H) 8q2u - 2q2 {l)\q q /2j ^(fn^J/lfn, fn) + Q 4 ^ plW^fJ 2 ■ 

We notice now that the middle contraction in the sum in (|4.1ip can only appear in the term p = q. 
Using the definition of G\f n and the second statement in (|4.10p we see that q 4 ql\\G q f n \\ 2 behaves 
asymptotically like 

^(moXa-i) 4 ' 1 ^" 1 ' 2 

Consequently, (|4.1ip has the same limit as 

(8<7 2 ^ - 4^ 2 ^!||^|| 2 ) -K 4^ 2 ^!||y„|| 2 - 2^ 2 (|) ! (^/ 2 ) ""^K^^^-^ ^*-> 



J Ufl_AlVf I' 1 Vll, T9/2, ,|2 



+ [q/2-l) H/n^ll^O, 

as n — > oo, where we have used the fact that ||/ n -^j^fnW 2 lj\ c q u i and {fn^J/^fn, fn) — > ^\ c q u - This 
proves the claim. □ 

We eventually deal with the convergence of the sequence 

MUM) = ~ f ^[{D z l {IqUn)> _ u} )D z I q (f n )\D z I q {f n )\} fi n (dz), n>l. 
Q J z 

Lemma 4.3. Under the conditions of Theorem \2.b\ we have that A^{I q {f n )) — > 0, as n — > oo. 

Proof. In view of the assumptions, we can directly apply Proposition 12.31 It follows that our claim is 
proved once we show that the three terms on the right-hand side of (|2.4p (with F = F n = I q (f n ) and 
H = n n ) converge to zero as n — > oo. Since the first term equals Ai(F n ), by virtue of the previous 
Lemma 14.11 we only have to prove the convergence of the remaining two summands. Our starting 
point is the following representation of the quantity (D Z2 D zl F n ) 2 = q 2 (q — l) 2 I 2 _2{fn(zi, Z2, • )), which 
is obtained by means of the product formula (|3.5p . Indeed, 

9-2 r / _ 9\ 2 / \ 
(D Z2 D zl F n ) 2 = q 2 (q-l) 2 Y J Y, rl ( q r ) ( p) h{ q -2)-r-ltfn{z u z 2 , ■ ) 4 f n (z u z 2 , •)) 

= q 2 (q-l) 2 I 2 (/„(*i, z 2 , •))• 
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Combining this representation with an iterated application of the triangle inequality, as well as of the 
isometric properties of multiple integrals, one deduces that the quantity 



^E JjD Z2 D Zl F n )* Mn (d*iK(dz 2 ) 
is bounded by a linear combination (with coefficients not depending on n) of quantities of the type 




||/n0l,*2,-) *r fn{zi,Z 2 ,-W ^n(dzi)l2 n (dz 2 ) = ||/n*„_j fn\\ ~> 0, 

IZJZ 

where the equality follows from a standard application of Fubini's theorem, and the convergence to 
zero is a consequence of the fact that a := q — £ 6 {2, . . . , q} and b := q — 2 — r 6 {0, . . . , a — 2}, as 
well as of the elementary identity \\f n *° / n || = ||/ n *g~ a / n || (2 < a < g). To deal with the remaining 
middle term, we use Fubini's theorem and the Cauchy-Schwarz inequality to deduce the estimate 



IE 



with 



J z Jz 

C n :=E £ (Jjp z2 D Zl F n f^ n {dz2)\ Uni&zi). 



Using again the explicit representation of (D Z2 D Zl F n ) 2 and applying several times Fubini's theorem, 
one sees that C n is indeed equal to a linear combination (with coefficients not depending on n) of 
objects of the type 

\\fn* b a /n|| 2 , with a = 2, . . . , q and 6 = 0,..., a -2. 
The conclusion follows immediately since our estimates do not involve the middle contraction. □ 

4.4 Proof of Proposition 12.91 

The product formula (|3.5p shows that 

(4.12) J 2 2 (/ n ) = / 4 (/n^/n)+4/ 3 (/n??/n)+/ 2 (4/ n ^/ n + 2/2) + 4/ 1 (/ n 4/ n )+2||/ n || 2 . 

Using the relation 

(4.13) 4!||/ n ^/ n || 2 = 2(2||/ n || 2 ) 2 + 16||/„ *J / n || 2 

(see e.g. formula (5.2.12)]), exploiting the orthogonality of multiple integrals of distinct orders 
and using the fact that ||/ 2 || — > by assumption, we infer that E[/|(/ n )] — 12E[/f(/n)] has the same 
limit as 

16 x 3!||/ n ??/ n || 2 + 16||/ n 4 f n \\ 2 + 48||/ n *{ f n \\ 2 - 96(/„ *} f n , f n ) + 3(2||/ n || 2 ) 2 
= 16 x 3!||/ n ^/ n || 2 + 16||/ n 4 f n \\ 2 + 48||/ n *} /„ - / n || 2 - 48||/ n || 2 + 3(2||/ n || 2 ) 2 . 

The conclusion follows by observing that ||/n|| 2 - > v by assumption, and then by applying Theorem 

Ym □ 
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4.5 Proof of Theorem I27T51 



Proof of Part A. According to Lemma [2. Ill since each F' n is completely degenerate, one has that F' n = 
h(fn), where /„ = h n /a(n), and the double integral is performed with respect to the compensated 
Poisson measure fj n = r] n — fj, n . It follows that the estimate (|2.12|) is a direct consequence of |22|, 
Theorem 4.2]. Using formulae (|4,12[) and (|4. 13j) . we deduce that 

E[7|(/ n )] = 16x3!||/ n ^/ n || 2 + 16||/ n 4/ n || 2 + 16||/„*l/ n || 2 
+ 2||4/ n ^/ n + 2/ 2 || 2 + 3(2||/ ri || 2 ) 2 , 

where the norms and contractions are of course taken with respect to the measure /%. Since 3(2||/ n || 2 ) 2 
converges to 3 by assumption, we deduce that, if (|2.1ip is verified, then the right-hand side of (|2.12p 
converges to zero, and therefore F' n converges in distribution to N. To conclude, observe that the 
estimates contained in [61 pp. 744-745] yield that E[(F n — F n ) 2 ] = 0(ra -1 / 2 ) as n — > oo, so that the 
estimate (|2.13p follows from the elementary inequality 

d w {F n , N) < d w (F n ,N) + [E(F n - Fn) 2 } 1 / 2 . 



Proof of Part B. Again in view of Lemma 12.111 and of the complete degeneracy of each F' n , we deduce 
that F' n = IzQin), where the double integral is again with respect to the compensated Poisson measure 
corresponding to rj n . The estimate (|2.15p is therefore a consequence of Theorem 12. 6\ and the fact that 
the distribution of F' n converges to T u is a direct consequence of Proposition 12.91 in the case h n = f n . 
The conclusion follows once again from the fact that E[(F n - F n ) 2 } = O (n x / 2 ) as n — > oo, in such a 
way that (|2.16p follows from the triangle inequality 

d 3 (F n ,T u ) < dz{F' n ,T u ) + [E(F; t - Fn) 2 } 1 / 2 . 

This completes the proof. □ 



4.6 Proof of Theorem I2TT71 

According to Theorem 7.3], one has that 

F'n = h{h n ) + Rn, 

where h n = h^ln = n^ 1 Ya=i e « ® e «> the double integral is realized with respect to the compensated 
Poisson measure fj n = rj n — nfi, and R n is a residual sequence of random variables such that 

E[R n } = 0(l/n), as n -> oo. 

It is immediate to verify that: (a) ( f z h n d/i 2 ) 1//4 = 0(1 /y/n) as n — > oo, (b) h n *\ h n = h n (where the 
contraction is realized with respect to fj, n ), (c) \\h n *\ h n \\ = 0(n~ 1 / 2 ) as n — > oo (since h n *\ h n (x) = 
n _1 Yli=i e i( x ) 2 )- The estimates are therefore a consequence of Theorem 12.131 (B). as well as of the 
estimates E[(F n — F n ) 2 } = 0(n -1 / 2 ) as n — > oo and 

d 3 (F n ,T u ) < dz(F' n ,T v ) + [E(F n - Fn?} 1 ' 2 - 

This completes the proof. □ 



4.7 Proof of Theorem I2T201 

We start with some general preliminaries which will be specialized below. Let F^ , . . . , Fn be centered 
square-integrable functionals of the Poisson measure r\ in the domain of the derivative operator D. 
For i G {1, . . . , d} let us define 



(4.14) 



a$ := E|2(F« + ^)+ " (DF®, -DL^F^ +E J ^F^D^F^^ 

+ E^( J D 2 l {F ( l)> _^ } )( J D,FW)| J D 2 L- 1 F«|^ n (dz), 
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and for i 7^ j G {1, ... , d} put 

(4.15) f}M := n{DF^,DL- x F^)l 7 W> := E J ^F^D^F^ Mn (ck). 

We estimate the distance between (the law of) F„ := (Fn\ . . . , i*n ) and (that of) T := (Gi, . . . , G^) 
by d(F n ,T) = sup |E0(F„) — E</>(r)|, where the supremum runs over all functions (/> : M d — >• IR whose 
partial derivatives up to order 3 are bounded, continuous and satisfy || • ||oo < 1. We notice that if 

d(F n , r) — > then F„ -4 T as n — > 00. 

Lemma 4.4. There exist constants K\ and K 2 such that 

d d 

d(F n , V) < K X °® +K 2 Y, [P ( n' j) + 7F) • 

i=l 

Proof. The technique adopted here is similar to the one used in the proof of the main result of [2] . To 
keep the argument more transparent and the formulas simpler we restrict ourselves to the case d = 2, 
the general case can be dealt with similarly. So, F n = (F^, F^) and T = {G\,G 2 ) and we have to 
show that 

(4.16) d((F« , if) ) , ( Gl , G 2 )) < K\ («« + a^)+K 2 +^ 2) W^ ] ) ■ 

To accomplish this task, we shall provide uniform estimates on \E(j)(F^\ F^) — Ecf)(Gi,G2)\- First 
write 

|E0(F«,F( 2 )) -EcP(G 1 ,G 2 )\ < |E[0(F«,F( 2 ))] -E^if ))]| 

+ \e[<p(g 1 ,fW)]-e[<p(g 1 ,g 2 )}\ = : |ri| + |T 2 |. 

We first deal with T 2 . Conditioning on G±, we are in a one-dimensional situation and can proceed as 

i<2\ 

in the proof of Theorem 12.11 This shows that T 2 contributes the term to the bound (|4.16p . We 
now consider the term T\ and write Cjj for the law of a random object U. Rewriting yields 



\6{x, y) - J (f)(g, y) C Gl (dg)^j ^ 2)) (d(x, y)) . 



For fixed y we consider the term in brackets as the left-hand side of a Stein-equation for the re- 
distribution so that 



U){x, y)- j (f){g, y) C Gl (dg)^j C (F m F (2 )) (d(x, y)) 

2(x + ui) + h' y (x) - xh y (x) C^ p (i) F ( 2 )^(d(x,y)) , 



(4.17) 



where, for fixed y, h y (x) is the solution of the Stein-equation associated with the test function x 1— > 
4>{x,y). We now consider the bivariate function h(x,y) := h y {x). Using the smoothness assumptions 
on (p together with the explicit representation 



h(x,y) 



-i(0(x,y)-E^(G(^),y)]) : x < -v x 

2{x+u 1 ) +gvi { x ) F- Vx {Hz,y) -E[</>(G(vi),y)\) gvi {z)dz : x > -v x , 



(recall the discussion preceding the proof of Theorem 12.11 and notice that g Ul ( • ) stands for the density 
of the law T Ul ) we deduce the following facts: (i) the mapping x \-t h(x,y) (for fixed y) is twice 
differentiable on v} (and it also admits right and left first derivatives at x = —v\ and (ii) 

the mapping y \-± h(x, y) (for fixed x) is twice differentiable on M.. All the involved derivatives are 
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bounded by a finite constant only depending on v\. Note that, in order to establish the estimates on 
y h-» h(x,y), one has to take derivatives under the integral and expectation signs, which is allowed 
thanks to the assumptions on <p. 

After these technical considerations we observe that (j4.17j) may be expressed in terms of h as 

2{x + vi) + dih(x,y) - xh(x, y) £ (F a)^ 2 ) ) (d(x, y)) 

(418) = E[2(^) + vi) + dMFk l \FV) - F x h{FV,FW)] 

= E[2(F« +v l ) + d 1 g{FP,F^) - {Dh^^-DL^F^)], 

where d\ stands for the partial derivative with respect to the first coordinate and where we have 
applied the integration by parts formula (|3.4[) of Malliavin calculus in exactly the same way as in the 
proof of Theorem 12. 1[ Using the notation Fn} z (v) = Frfi (v + S z ) — F^ (rj) for i 6 {1, 2} and z G Z, we 
may write 

D z h(F 1 ,F 2 ) = h(F(%F$) - h{F£\FW) 

= (kFl%F$) - h(F(%FP)) + (h(Fi%FP) - h(FW,FW)) =: S t + S 2 

Thanks to the properties of h described above, we find that 

Sx = d 2 h(FS,F^)D z F^ + rW(D m fW>) and S 2 = d 1 h{F^ y F^)D z F^ + R^(D zF (% 

where and are such that 

\R^{D Z F^)\<K^\D Z F^ 

and 

where d\\ and d 22 , respectively, denote the second derivative with respect to the first and second 
coordinate and where K j 1 ^ , K , iiQ are finite constants. Combining this with (|4.18p and taking the 
supremum over all <fi, we obtain the contributions a„ , fin' 1 " 1 and jn' 1 ^ in (|4.16p . Inverting the role 
of F^ and Fn" 1 in the previous discussion gives the bound (|4.16p . with constants K\ and K 2 only 
depending on {i>\, v 2 ). □ 

Proof of Theorem \2.20L Let us define the random vector I n := [l qi (fn ),•••, ^ (/n )) • We shall prove 
that d(I n ,r) — )• as n — > oo. Lemma 14.41 implies that for this it is sufficient to check that 0$ — > 0, 
/3 n iJ) -»■ and that 7^ -> as n — >• 00 for any combination of i and j. Under the assumptions 
in the statement, writing Fn = Iq^fn ) one has the following three facts for every £ = 1, . . . ,d: (a) 

a® — )• 0, as n — > 00, (b) as n — > 00, 

E J (D Z F®) 4 f, n (dz) -+0, 

and (c) the sequence 

E / (D Z F^) 2 /j n (dz) = g?E / (DL~ l F^) 2 fi n (dz), n > 1, 

is bounded. An application of the Cauchy-Schwarz inequality yields therefore that 7^ — > for any 
allowed choice of i and j. To check the fact that j3 n —> 0, we apply once more the Cauchy-Schwarz 
inequality to obtain 

2\ V2 



/?F < 9i 2 (e^I^ ■)) -))^(dz 
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We use now the general product formula (|3.5p for multiple integrals to express 

I«-i(fk%,'))l«-i(fP(*,')) 

as a sum of multiple integrals and the stochastic version of Fubini's theorem allowing us to exchange 
deterministic with stochastic integration; see [23} Theorem 5.13.1]. By assumption, < qj. Using the 
triangle inequality several times yields 



with the constant K(r, £, qi, qj) given by 

Kirx^) =(r- ±v( q ;:i) ( q i:l) (; : ;) c* + * - r - ^ 

The proof is completed by observing that (see [241 Lemma 2.9]) 

II f(i)+t f(J)\\ < II f (0 J- f(*)|| || fU) J- f(i)|| v n 

WJn *rJn II — HJn Jn II IUn *r Jn II ^ u 

for all choices of i, j, because of the assumptions in the theorem and the fact that 2qi ^ qj for i ^ j. □ 
4.8 Proof of Theorem 

We start again with some preliminaries. Let F^, ■ ■ ■ , Fn dl+d2 ' > be square integrable functionals of the 
Poisson measure r\. For i G {di + 1, . . . , d\ + cfe} let us define 

<5« :=E|1- {DFi i \-DL- 1 F^)\+¥. f JD Z F^\ 2 \D Z L- 1 F^\ ^ n {dz) 

and for i G { 1 , . . . , <ii } let a['' be as in (|4.14p and for i,j G {1, + ofo} let and 7n be 

as in (|4.15p . We will estimate the distance between (the law of) F n := (F^\ . . . , Fn dl+d2 ^ and 
(that of) the hybrid vector H := [Gi, . . . , , N^+i, ■ ■ ■ , Nd 2 ) by the hybrid distance dh(F n , H) = 
sup|E(/>(F n ) — E0(H)|, where the supremum runs over all functions <p : M. dl+d ' 2 — > R whose partial 
derivatives up to order 3 are bounded, continuous and satisfy || • ||oo < 1. 

Lemma 4.5. There exist constants K\, K% and K3 such that 

d h (F n ,H) <K,J2 ^ +K 2 6 n + K ^H (P ( n' j) + • 

i=l i=d±+l i,j=i 

Proof. This follows along the same lines of argumentation as the proof of Lemma 14.41 For this reason 
the details are omitted. □ 

Proof of Theorem \2.21\ We first use Lemma 14.21 to see that because of (|2.17p . 0$ — > for any 
i G {l,...,di}. Next, we apply [221 Theorem 5.1] to infer that under (127T8J) . 5^ -> as n — > 00 
for any i G {d±, +1, . . . ,d\ + c^}. The remaining discussion of fin'^ and 7n is very similar to the 
multivariate pure Gamma case so that -)• and 7^ j) -> for all* ^ j G {1, . . . , d x + d 2 }. In 

view of Lemma 14.51 this completes the proof. □ 
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4.9 Proof of Theorem 12^51 

We consider a measurable bounded test function <j) : R x Z + — > R such that <j> has uniformly bounded 
derivatives up to the order three in the first variable. By a slight variation of the arguments leading 
to the proof of [21 Theorem 2.1] one has that there exists a universal constant K > (independent of 
n) such that 

|E[0(/ g (/ n ), H n )) - E[<t>(G, P)]\<K (A n + B n + C n + D n ), 
where (similar to a n etc. above) 

A n := E \2(F n + u)+ - (DF n , -DL~ x F n )\ + £ E[\D z F n \ 2 \D z L- x F n \] fi n (dz), 

B n := \E[H n ] - A| + E | A - {DH n , -Dir x H n ) \ + ^ E\D z H n (D z H n - \)D z L- x H n \ » n (dz), 
C n :=E[(\DH n \,\DI q (f n )\)} 

and 

D n :=E^( J D 2 l {/9(/n)> _, } )(D z / (? (/ n ))|^L- 1 F n |^ n (dz) 

In view of Theorem 12.61 (as well as of the estimates leading to its proof), the assumptions in the 
statement imply that A n + B n + D n — > 0, and, moreover, that 



E J (D z I q (f n )) 4: fj, n (dz) — > as n -> oo. 



The conclusion is obtained by observing that, by virtue of Holder's inequality, and since DH n takes 
values in Z, 

C n < (e jjD z I g (f n )) 4 ^(dz)^ ' x (e j^D z H n ) 4 / 3 fj, n (dz)\ ' 
< (e J^(D z I q (f n )) 4 fi n (dz^j ' x (e J^(D z H n ) 2 ^(dz)^ 1 
where we have implicitly used assumption (|2.19p . □ 
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